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ABSTRACT 

We present results from a deep photometric study of the rich galaxy cluster Abell 
2218 (z = 0.18) based on archival HST WFPC2 F606W images. These have been 
used to derive the luminosity function to extremely faint limits (A/f606W ~ —13.2 
mag, no « 24.7 mag arcsec -2 ) over a wide field of view (1.3 h~ 2 Mpc 2 ). We find the 
faint-end slope of the luminosity function to vary with environment within the cluster, 
going from a = —1.23 ± 0.13 within the projected central core of the cluster (100 < 
r < 300 hr 1 kpc) to a = -1.49 ± 0.06 outside this radius (300 < r < 750 hr 1 kpc). We 
infer that the core is 'dwarf depleted', and further quantify this by studying the ratio 
of 'dwarf to 'giant' galaxies and its dependency as a function of cluster-centric radius 
and local galaxy density. We find that this ratio varies strongly with both quantities, 
and that the dwarf galaxy population in A2218 has a more extended distribution than 
the giant galaxy population. 
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The galaxy luminosity function (LF) - the number density 
of galaxies per unit luminosity interval - is a fundamental 
descriptor of the galaxy population and, as such, contains 
important information on the formation and evolution of 
galaxies (Bingelli, Sandage & Tammann 1988, Benson et al. 
2003). At low redshift, the LF will contain the combined im- 
prints of the galaxy initial mass function (Press & Schechter 
1974, Schechter 1976), together with the effects of any sub- 
sequent evolutionary processes which modified this distri- 
bution (e.g., hierarchical merging; White & Rees 1978). De- 
terminations of the LF in different environments and at a 
variety of redshifts provide the only direct hope of disentan- 
gling these environmental and evolutionary effects. 

In this context, rich clusters are fundamental testing 
grounds, representing the densest environments in which 
galaxies reside and, within the hierarchical clustering frame- 
work, the ultimate examples of where build-up through 
merging has taken place. Here galaxies at all luminosities, 
both bright and faint, may bear evidence of environmen- 
tal effects: the 'giant' (L > L*) galaxies at the bright end 
of the LF are the strongest candidates for having been 
built up through successive mergers and accretion (e.g. 
Kauffmann, White & Guiderdoni 1993). Equally, the sub- 
luminous (L < L*) 'dwarf galaxy population may be a vic- 
tim of, and therefore depleted by, these merging and can- 



nibalisation processes. In contrast, galaxy 'harassment' may 
boost the dwarf galaxy population through whittling down 
more luminous cluster galaxies as they undergo high speed 
encounters (Moore et al. 1996). Additionally, both giant and 
dwarf galaxies in the core regions of clusters are likely to 
be susceptible to such processes as ram-pressure stripping 
(Gunn & Gott 1972) and cluster tidal effects (Byrd & Val- 
tonen 1990, Bekki et al. 2001), which severely modify their 
star forming activity and hence their luminosity (see also 
Bower & Balogh 2003, Moore 2003). 

This potential for the rich cluster LF to be used as an 
evolutionary 'probe' has led to numerous efforts to measure 
it. In recent times, the faint-end slope of the LF - described 
by the parameter, a, in the Schechter (1976) function rep- 
resentation - has been of particular focus, with numerous 
claims that it is significantly steeper in clusters than in the 
general field (De Propris & Pritchet 1998, Trentham & Tally 
2002). 

Comparison of cluster LF results is, however, compli- 
cated by a number of factors: the differing magnitude ranges 
over which the slope is calculated, the different methods of 
field galaxy subtraction employed by different authors, and 
the different physical cluster cross-sections covered (Driver 
& De Propris 2003). A further complication is that some 
observations suggest the numbers of dwarf galaxies varies 
with environment - what Phillipps et al. (1998) referred 
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to as a "dwarf population-density" relation. This manifests 
itself as the densest regions in clusters having the lowest 
dwarf galaxy fractions, the ratio of dwarf to giant galax- 
ies increasing with decreasing giant galaxy density, and the 
dwarf populations in clusters being more spatially extended 
than the giants (Smith et al. 1997, Driver et al. 1998a, 
Driver et al. 2003). This is also compounded by the like- 
lihood that there exist two dwarf populations, such as seen 
in Virgo, whereby the dwarf ellipticals are centrally concen- 
trated (e.g., Binggeli, Sandage & Tammann 1985; Conselice 
et al. 2003) and the dwarf irregulars (and/or dwarf low sur- 
face brightness galaxies) lie predominantly in the halo (Saba- 
tini et al. 2003). 

In this paper we present a case study of the well 
known rich cluster Abell (A)2218 at z = 0.18, deriving 
its LF to very faint limits and over a broad range of en- 
vironment from deep, wide-field imaging obtained with the 
Hubble Space Telescope (HST). Famous for its spectacular 
strong gravitational lensing of background galaxies (Kneib 
et al. 1996), A2218 is an Abell richness class 4 cluster 
with a Bautz Morgan type of II, whose centre is domi- 
nated by a large, low surface brightness cD galaxy with 
an envelope extending over more than 180h _1 kpc (Pello- 
Descayre et al. 1988). It also has a massive, deep potential 
well, with a central velocity dispersion of 1370 km s -1 (Le 
Borgne et al. 1992), and is a strong X-ray emitter [Lx(0.5- 
4.4keV)= 3.3 x 10 44 h~% ergs^ 1 ; Jones et al. 1993]. More- 
over, a strong Sunyaev-Zeldovich decrement is observed in 
its direction (Birkinshaw & Hughes 1994). As well as pro- 
viding an extreme case in the context of observing envi- 
ronmental effects, A2218 is also a propitious choice in that 
its high richness makes it particularly amenable to reliable 
'background' removal using statistical methods (Driver et 
al. 1998b). 

The layout of this paper is as follows: In section 2 
we describe the HST imaging of A2218, and the proce- 
dures followed in the detection and photometry of objects, 
star/galaxy separation, and the extraction of the cluster 
population through the careful subtraction of the back- 
ground population. For the latter, considerable attention is 
paid to the numerous effects that contribute uncertainties to 
this process, in particular the effect of gravitational lensing. 
In section 3 we present the luminosity function, subdivided 
into inner (projected core) and outer (halo) components. We 
then determine the dwarf-to-giant ratio versus giant galaxy 
density and investigate the radial profiles of the giant, dwarf 
and ultra-dwarf galaxies (defined below) . We summarise and 
discuss our results in section 4. Throughout we adopt a 
Qm = 0.3, n A = 0.7 and H = lOOkms^Mpc -1 cosmol- 
ogy; this puts A2218 at a distance modulus of 39.01 mag 
(inclusive of a fc-correction), with 1 arcmin projecting to 128 
co-moving kpc at this distance. 

2 DATA ANALYSIS 

Our study is based on archival HST 1 Wide Field Planetary 
Camera 2 (WFPC2) images of A2218, obtained as part of 

1 Based on observations made with the NASA/ESA Hubble 
Space Telescope, obtained from the data archive at the Space 
Telescope Science Institute. STScI is operated by the Associa- 
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Figure 1. The HST WFPC2 mosaic used for this study, showing 
the arrangement and location of the 22 pointings on the sky. A 
total exposure time of 8,400 s was obtained at each pointing. This 
is overlayed on top of a digital sky survey image of the A2218 field. 

a program to study weak gravitational lensing (PI: Squires, 
PID:7343) over a much larger region of this cluster than 
that covered in the original strong-lensing study (Kneib et 
al. 1995). These data comprise a mosaic of 22 pointings in 
the F606W passband. Each pointing involved a total expo- 
sure time of 8,400 s. The location of these pointings on the 
sky is shown in Figure 1; the mosaic provides coverage of 
a roughly square field with a total area of ~ 82 arcmin 2 . 
The resolution, depth and coverage of the data are unique 
and particularly well suited to investigating A2218's galaxy 
population down to very faint limits and out to a radius of 
~ 0.75 h~ 1 Mpc from its centre. 

The HST data were obtained from the HST archive 
with the standard procedures of bias removal and flat-field 
division having been applied. Each pointing consists of 12 
dithered exposures, each of 700 s duration. In order to re- 
move cosmic rays, the individual frames were drizzled to- 
gether using the drizzle task in the IRAF dither package, 
following the procedure of Koekemoer et al. (2002) . The fi- 
nal drizzled image had a pixel size of 0.0498 arcsec. Only 
the portions of each pointing with the full 8,400 s of expo- 
sure were used in the analysis, resulting in a final overall 
field-of-view of 0.0228 sq. degrees or 1.34 Mpc 2 at the clus- 
ter redshift. The photometric zero point adopted for our 
data was that published for the WFPC2 by Holtzman et al. 
(1995), which places the instrumental F606W magnitudes 
onto the VEGA system. The quoted accuracy of this zero- 
point is 2%. 



2.1 Object detection and photometry 

Objects were detected and photometered using the SEx- 
TRACTOR package of Bertin & Arnouts (1996). Based on 
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Figure 2. A typical image obtained with one of the WFPC2 
CCDs (WF3) with SExTRACTOR detections circled. The image is 
80" (0.17 h,- 1 Mpc) on a side. 



previous experience in using SExTRACTOR with WFPC2 im- 
ages, a detection criterion of 20 connected pixels 1.5a above 
the sky RMS was used. Figure 2 shows a typical image 
obtained with one of the WFPC2 CCDs (WF3), with the 
objects detected by SExTRACTOR circled. The magnitude 
best_mag was used as our measure of total galaxy magni- 
tude (hereafter denoted simply as F606W). This adopts a 
Kron magnitude as the default value, except for crowded 
regions where an extrapolated isophotal magnitude is used. 
Note that radius of the Kron (1980) extraction aperture is 
set to 2.5 times Rk where Rk is the first moment of the im- 
age distribution. A total of 8038 objects were detected with 
F606W < 25.8 mag (equivalent to Mf606w < —13.2 mag), 
the magnitude at which the number of detections turned 
over. 

The SExTRACTOR parameter CLASS_STAR was used as 
the basis for separating stars from galaxies, supplemented 
by visual inspection in the small number of cases where the 
software was equivocal in its classification. 

A broad overview of the range and depth of the ob- 
jects detected by SExTRACTOR and in particular the pop- 
ulation of galaxies included in this study, is shown in Fig- 
ure 3. Here we plot the mean central surface brightness, 
/iapp (measured in a circular aperture of area equivalent to 
the 20-pixel minimum required for object detection), of all 
our detected objects as a function of their apparent total 
magnitude. This reveals three distinct populations of ob- 
ject: cosmic rays, which form the tight upper-most diago- 
nal locus, stars, which form the slightly broader central lo- 
cus, and galaxies, which form the lower-most broad wedge 
of objects. The stars and galaxies are clearly distinguish- 
able down to F606I4 7 ~ 24. Also of note is the surface- 
brightness 'edge' to the galaxy population, with no objects 
seen fainter than /x apP ~ 24.5 mag arcsec 2 . This is consistent 
with the 1.5cr isophotal detection limit used in SExtrac- 
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Figure 3. The aperture surface brightness of all our detected ob- 
jects, plotted as a function of apparent magnitude. The measured 
objects, represented by the dots, are clearly seen to define three 
distinct loci within this diagram: cosmic-rays (upper- most locus), 
stars (central locus), and galaxies (lower wedge). Stars and galax- 
ies are easily separated down to a magnitude of F606W ~ 24.5. 



TOR, which corresponds to fi « 24.7 mag arcsec -2 . Given 
the way the galaxy wedge broadens with increasing magni- 
tude and surface brightness, it is clear that our sampling of 
the galaxy population is surface-brightness limited at mag- 
nitudes fainter than F606W ~ 23. However, as we discuss in 
the following section, this has little impact on our ability to 
accurately quantify the cluster galaxy counts at these faint 
apparent magnitudes. 



2.2 Field galaxy removal 

In the absence of spectroscopy, the removal of the foreground 
and background 'field' galaxy populations must be done sta- 
tistically. Here we follow the strategy adopted by Driver et 
al. (2003) for a similar HST-based study of the rich clus- 
ter Abell 868 (z = 0.15). This involved estimating the field 
galaxy counts from the deepest available HST/WFPC2 data 
taken in the F606W band. This comprised images of the 
Hubble Deep Field (HDF) North (Williams et al. 1996), the 
HDF-South (Casertano et al. 2000), and 10 other fields with 
total exposure times between 10,000-20,000 s (and hence of 
comparable depth to our A2218 cluster mosaic). Object de- 
tection and photometry was performed on these images us- 
ing SExTRACTOR, with an identical set of parameters to 
those used on our A2218 data. 

The galaxy number counts derived from these data are 
shown in Figure 4, where we have plotted the combined 
counts from the two HDF fields and those from the 10 other 
reference fields separately. To within the uncertainties, we 
see that there is no significant difference in the galaxy counts 
between the 'deep' (HDF) and 'shallow' reference fields. This 
is important, since by including the much deeper HDF data 
(surface brightness limit of p apP ~ 26.0 mag arcsec 2 cf. limit 
of At app ~ 24.5 mag arcsec for our 'shallow' reference fields 
and A2218 - see Fig. 3) in our field galaxy count estimates, 
we want to be certain that this does not lead to an over- 
estimation through the inclusion of an additional population 
of faint, low surface brightness galaxies unseen in our cluster 
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Figure 4. The combined HDF-North and HDF-South galaxy 
counts (open squares) and the combined shallower reference 
counts (closed circles), plotted on a linear scale to aid comparison. 
The consistency of the counts implies negligible selection bias in 
either the detection or photometry between the shallow and deep 
fields. 




Figure 5. Top panel Galaxy number counts as a function of 
F606W magnitude for Abell 2218 (diamonds), the combined 
HDF-North, HDF-South and WFPC2 counts (open circles) and 
the MGC counts (triangles). The solid line is a polynomial fit to 
all the reference field counts. Bottom panel The galaxy number 
counts shown relative to the N <x lo°' 4m relation. 



images. The agreement seen between the two different sets 
of number counts shown in Figure 4 would indicate this is 
not the case. Hence we can be confident that our selection 
limits are identical both inside and outside the cluster, and 
the excess galaxy counts we observe after field subtraction 
is due to having detected a bona fide population of cluster 
galaxies rather than one that is spurious. 

In the top panel of Figure 5 we plot in the tradi- 
tional way, the galaxy number counts for both the 'field' 
(as per Figure 4) and what we measure in the direction 
of A2218. We also supplement the counts at bright magni- 
tudes (F606W < 17.5) by plotting those derived from the 
ground-based Millennium Galaxy Catalog (MGC) by Liske 
et al. (2003). We caution that this has involved transform- 
ing the MGC counts from the B-band (in which they were 
measured) onto our _F606VF(VEGA) system. This we did 
using the relation Bmgc - F606VF(VEGA) = 1.06 derived 
by Driver et al 2003. This simple colour conversion is only 
strictly valid for the very bright end of the MGC counts: 
F606W < 17.5 mag. However, since their only purpose in 
this study is to provide a bright magnitude 'anchor' point - 
the accuracy of which has a negligible effect on the overall 
accuracy of our field galaxy subtraction at the much more 
critical fainter magnitudes - this is quite satisfactory. 

Comparison of the A2218 and field counts in the upper 
panel of Figure 5 shows a clear divergence between the two 
with decreasing magnitude, with the field counts having a 
slope close to the canonical value of ~ 0.4 observed in other 
studies at these wavelengths (e.g., Metcalfe et al. 2001). To 
provide a clearer picture of the deviation between the A2218 
and field counts, we plot in the lower panel of Figure 5 the 
counts relative to the N oc io 0,4m relation (where m repre- 
sents our F606W magnitudes), appropriately normalised to 
the A2218 mosaic field of view. This reveals some residual 
magnitude dependence in the field counts, so to account for 
this and ensure that they are accurately represented over 
the entire magnitude range covered by our study, we fit a 
quadratic to the combined MGC+WFPC2 counts, yielding: 

dl °f loiV = -12.9346 + 0.96077m - 0.012874m 2 . (1) 
dm 



This quadratic fit is shown as the solid line in both panels 
of Figure 5. 

Table 1 provides a tabulation of the galaxy counts (and 
their associated errors) relevant to our A2218 analysis. In 
column 2 we list the total counts, n tot , observed within our 
A2218 field and in column 5 our estimates of the numbers of 
field galaxies, n re f, within the A2218 mosaic field-of-view, 
based on the above quadratic fit. Columns 3 and 6 give esti- 
mates for the Poisson noise in the components that make up 
each observed count: the standard deviation of the counts 
along the cluster line of sight: 

ftot = V n tot, (2) 

and the error due to Poisson noise in our estimate of the 
mean field (or reference) count: 

n / ^A2218 Href /q\ 

^ rcf _ V «wppc2 12 ' [6) 

respectively. Here Oa2218 and £>wfpc2 represent the area 
of sky covered by our A2218 observations and within one 
WFPC2 field, respectively. 

There is an additional uncertainty in the galaxy counts 
as a result of general galaxy clustering along the line of sight 
to A2218 (Glazebrook et al. 1994, Djorgovski et al. 1995). 
The error introduced by this effect has been computed by 
Huang et al. (1997) and is given by their equations 5 and 6. 
An analogous formula is also given by Driver et al. (2003; 
equation 7) . This is a vital component which typically dom- 
inates the error budget at faint magnitudes and has often 
been neglected in previous studies. The errors introduced 
into our counts by general galaxy clustering along the line 
of sight to the A2218 field (C c ) and to the reference fields 
(C rc f) are given in columns 4 and 7, respectively. They have 
been calculated using the Driver et al. equation: 

e *(N) = (^)-°- 8 ■ n- aA ■ N 2 (m) • i(r - 235m + 2 - 73 . (4) 

Here, N is the field number counts (per 0.5 mag interval) 
and fi is the area of sky over which the counts are mea- 
sured. The field galaxy counts appropriate to our A2218 
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Figure 6. Comparison of the measured variance in the 12 refer- 
ence fields (filled circles) and the predicted variance as described 
in the text (solid line). The dashed line and dotted line show the 
Poisson and clustering contributions to the predicted variance re- 
spectively. Note that the error bars plotted for the observed points 
were determined using a 'jack-knife' method. 



field (and hence to the derivation of the 'clustering' er- 
rors listed in column 4) are given in column 5 of Table 1. 
These have been derived using the above quadratic repre- 
sentation of the MGC+WFPC2 counts, with the appropri- 
ate normalisation applied to account for the different sky 
coverage: Oa2218 = 0.0228 sq. deg, Omgc = 30.310 sq. deg, 
ftwFPC2 = 0.0011 sq. deg. The small area of the WFPC2 
field means that the counts observed within it are more 
susceptible to line-of-sight clustering, as can be seen in the 
larger 'clustering' errors listed for the field in column 7 of 
Table 1. 

In Figure 6 we show how the observed variance in the 
galaxy number counts from field to field compares with what 
we would expect based on our above error budget analysis. 
It can be seen that, overall, there is good agreement between 
the two, with our predictions, if anything, being slightly pes- 
simistic. 

The final column of Table 1 contains the field- 
subtracted 'cluster' counts within our A2218 field, with their 
quoted uncertainties representing the combination of all the 
contributing errors described above. Figure 7 shows graphi- 
cally the contribution of each of the four errors to the final 
error budget. 



2.3 Lensing 

A further possible correction to the background counts con- 
cerns the effects of gravitational lensing of the background 
population by the cluster mass (Trentham 1998). This con- 
sists of two competing effects: a magnification effect, where 
the background population is brightened and therefore seen 
in greater numbers at any given magnitude, and an apparent 
expansion of the background field behind the cluster, which 
reduces the surface density of background objects. The pre- 
cise details of the combined correction required to account 
for these two effects are given in Bernstein et al. (1995) and 
Trentham (1998); we therefore give here only a brief out- 
line of the model which we will use to estimate the effect of 
lensing in A2218. 




Figure 7. The errors in the Abell 2218 number counts as a 
function of magnitude. The four individual error components are 
shown along with the combined total error (solid line). 



The quantity that requires computation is the fractional 
change in the counts due to lensing: 



N'(m) 
N(m) 



(5) 



where N(m) is the observed reference counts per unit area 
and N' (m) is the counts we should instead observe when the 
lensing effect of the cluster is taken into account. Calculation 
of N'(m) requires knowledge of the number of galaxies per 
unit redshift per magnitude per steradian, which is given by: 



. ln(10) dV 

n(m,z) = — — 4>(m,z) — , 
llm az 



(6) 



where 4>(m, z) is the field galaxy luminosity function at red- 
shift z. We assume that this can be described by a Schechter 
function at all redshifts, with luminosity evolution parame- 
terized by: 

<f>*(z) = aS*(0)(l + zf 



M*(z) = M*(0) + 51og(l + z) 



(7) 
(8) 



(see Broadhurst et al. 1995, Trentham 1998). We use the lo- 
cal luminosity function measured from the 2dF Galaxy Red- 
shift Survey (Norberg et al. 2001) and the colour transfor- 
mations from Driver et al. (2003) which, after conversion to 
our cosmology, results in <f>* = 0.0169Mpc~ 3 , M* = -20.65 
and a = —1.21. The volume element is given by: 

— -nx 2 — (9) 

dz ~ llX dz> (9j 
where x is the comoving distance, and for a flat cosmology: 
dx c 
dz ~ H»yj£l m (l + zY + Q.A 

Having defined the model n(m, z) distribution, the true 
N(m,z) can be recovered by requiring that J Q N(m,z)dz 
reproduces the observed number counts, that is: 



(10) 



N(m,z) = 



N(m) 



-n(m, z) 



J o °°n(m, z)dz 
The function N'(m) is then given by: 

N'(m,z) = ^N(m + 2.5logA,z), 



(11) 



(12) 
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Table 1. Summary of the number counts for A2218. The total galaxy counts detected in the A2218 field 
(column 2) and the field reference counts (column 5) along with all the associated errors (see text). The final 
column gives the field subtracted galaxy counts along with the combined error. 



where A is the amplification of the lens. Again, following 
Trentham, we assume the cluster mass profile to be a spher- 
ically symmetric isothermal sphere. A is then given by: 



A(r,z) 



4trj 2 D c D C s(z) 



rc 2 D s (z) 



(13) 



where D c , D s and D cs are the angular diameter distances 
to the cluster, to the source and from the cluster to the 
source respectively (Miralda-Escude 1991). fi ens can then 
be calculated by radially averaging ^- over the annulus for 
which we want to evaluate the lensing correction: 



2 C /„" N(m, z)dz + /~ N' (m, z, r)dz 



rdr 



-,(14) 



flens(m) = 



where ri and V2 are the inner and outer radii of the annulus. 

In Figure 8 we plot /i cns for a series of annuli in which we 
will derive cluster LFs in our subsequent analysis. The solid 
line shows /i cns for the inner region of the cluster bounded 
by n — 100 kpc and = 300 kpc. It is this region where the 
effects of lensing are at their maximum. Note that the inner 
100 kpc has been excluded due to the unrealistic nature of 
the isothermal model in that region (AbdelSalam, Saha & 
Williams 1998). The dotted line shows /i ens for the outer 
region of the cluster, with boundaries r\ — 300kpc and = 
750kpc. Finally, the dashed line shows /i en s for an annulus 
encompassing the 'whole' field: r\ — lOOkpc to ri = 750kpc. 

In general, the family of curves shown in Figure 8 are 
in excellent agreement with those derived by Trentham (for 
Abell 665 at z — 0.18; see his Figure 3) in terms of their 
overall behavior with apparent magnitude, the amplitude 
of the effect (<; 10%), and its much diminished importance 
at larger clustercentric radii. The only systematic differ- 
ence of note is that our /i cns function drops below unity at 
faint magnitudes (F606W > 23) - that is, lensing decreases 
(rather than boosts) the number of faint background galax- 
ies seen in the direction of the cluster. This difference can 
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Figure 8. The function /i en s, averaged over three different annuli: 
lOOkpc < r < 300kpc (solid line), lOOkpc < r < 750kpc (dashed 
line) and 300kpc < r < 750kpc (dotted line). 



be attributed to slightly different assumptions made about 
the n(m, z) distribution, in particular the differences in the 
representations of the observed galaxy counts, N(m), in nor- 
malizing n(m,z). Here Trentham uses a linear function to 
represent d log N '/dm, whereas we have used a quadratic rep- 
resentation (see Figure 5). 



3 RESULTS 

3.1 Luminosity Function 

As a first step in our analysis, we take the final field- 
corrected cluster counts listed in the final column of Table 1 
and construct the luminosity distribution (LD) for A2218's 
galaxy population; this is shown in Figure 9 (filled circles) . 
Formally, this distribution is well described by a Schechter 
function with Mp 606W = -20.52±0.27 and a = -1.38±0.05 



over the range —21.5 < Mf606W < —13.2 (solid line in left 
panel of Figure 9). Because the errors in M* and a are 
tightly correlated, we present alongside our derived LD the 
error contours for these two parameters (the thick lines in 
the right-hand panel). Also shown in this figure is the LD 
recovered after the reference field counts have been scaled 
by /lens (dashed line in Figure 8) to correct for lensing effects 
(open circles). It can be seen that this makes very little dif- 
ference to our observed LD. This is reflected quantitatively 
in the values of the fitted Schechter function parameters, 
-Mpeoew = -20.25±0.28 and a = -1.40±0.05, which differ 
from the original values by less than la. 

For comparison, we also plot in Figure 9 the Mp 606W 
and a values for the Schechter function fits to the compos- 
ite LF derived for the z ~ 0.1 clusters in the 2dF Galaxy 
Redshift Survey (2dFGRS; De Propris et al. 2003), and the 
LF derived for the rich cluster Abell 868 (z = 0.15) from a 
comparable HST-based dataset to what we have used here, 
by Driver et al. (2003). It can be seen that our A2218 LF 
has a steeper faint-end slope (more negative value of a) and 
a slightly fainter M* than these other two LFs, particularly 
for the 2dFGRS clusters where, given the size of the er- 
ror bars, the differences are reasonably significant. On the 
other hand, the 2dFGRS and A868 LFs do not extend as 
faint as ours (Mf606W — — 16 cf. Mf606W — —13), and it is 
well known that the value of the parameter a is quite sen- 
sitive to the luminosity range covered, with the faint end 
steeping at fainter magnitudes (Lobo et al. 1997). A bet- 
ter comparison is obtained by fitting the A2218 data over 
the same absolute magnitude range, for which we recover 
M F606W = -20.26 ±0.27 and a = -1.23 ±0.09. This is also 
the range over which our sample is complete with respect to 
surface brightness. The error contours for the LF parameters 
fit over this range are shown as the thin contours in Figure 
9. We see that the LF in this case has a shallower slope, with 
its a value being consistent with those of the 2dFGRS and 
A868 LFs. However, the small offset in M* still remains. 
Such small differences may well result from the regions of 
the clusters being compared here not being identical; the 
likelihood of this being the case is highlighted below. 

Finally, it is also worth noting that we see evidence of 
an 'inflexion' in A2218's LF at Mf606W ~ —18, something 
which has been seen in many other clusters (e.g. Driver et 
al. 1998a; Barkhouse, Yee & Lopez-Cruz 2002). 

The wide field provided by the 22 WFPC2 pointings 
allows us to sub-divide the cluster into projected 'core' and 
'halo' components and re-derive the luminosity distributions 
for each component. To do this we split the inner and outer 
regions at 300/i _1 kpc from the central, dominant cD galaxy. 
Somewhat arbitrarily, this radius was chosen to be twice 
the typical cluster core radius, as measured, for example, 
for the EN ACS sample (Adami et al. 1998). In assembling 
the luminosity distribution for the inner region, we exclude 
the central lOO/i -1 kpc because of our inability to compute 
lensing corrections in this region (see Section 2.3). The outer 
radial limit is « 750/i~ 1 kpc. 

Figure 10 shows the resulting LDs with both the uncor- 
rected (solid symbols) and lensing corrected (open symbols) 
data points plotted. Our Schechter function fits and error- 
contours are shown in the usual way; because the lensing 
corrections have a negligible effect, only those based on the 
uncorrected data are displayed. From Figure 10 we see that 
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Figure 9. Left panel: The galaxy luminosity distribution derived 
for A2218. The data are represented by circles; those that are 
filled have had no correction for the effects of gravitational lens- 
ing, whereas the open circles include this correction. The data 
point contributed by the central cD galaxy is shown as the filled 
square. A Schechter function fit to the (uncorrected) data points 
is represented by the solid line. Right panel: The la, 2cr and 3a 
error contours for the Schechter function parameters M* and a 
for M < —13.2 (thick contours) and M < —16 (thin contours) . 
For comparison, the values of these parameters for the 2dFGRS 
composite cluster LF (Dc Propris et al. 2003) and that of A868 
(Driver et al. 2003) are represented by the open and filled circles, 
respectively. 

the LD for the central region has a significantly flatter faint- 
end slope than that for the outer region. The inner 'core' LF 
has a slope of a — —1.23 ± 0.13 while the outer 'halo' LF 
has a = —1.49 ± 0.06. This is in line with recent claims 
that the LF of galaxies in the Coma cluster is flatter in the 
vicinity of the two central dominant giants, with it steep- 
ening as a function of clustercentric distance (Beijersbergen 
et al. 2002; see also Lobo et al. 1997). Driver et al. (2003) 
obtained a a similar result from their LF analysis of A868. 
Barkhouse et al. (2002) also find a significant steepening of 
the LF as a function of clustercentric radius in their sample 
of local clusters as do Sabatini et al. (2003) for the Virgo 
cluster and Andreon (2002) for the z = 0.31 cluster A2744. 
Similarly, luminosity segregation may be responsible for the 
dwarf population density relation discussed in the following 
two sections. However, Mobasher et al. (2003) argue against 
the claim of Beijersbergen et al. (2002) for luminosity seg- 
regation in Coma, while Paolillo et al. (2001) fail to find 
steeper LFs in their composite LF for the outer regions of 
nearby clusters. 



3.2 Dwarf-to-giant ratio 

A simpler, non-parametric method to test for the luminos- 
ity segregation indicated in Figure 10, is to measure the 
dwarf-to-giant ratio (DGR - Driver et al. 1998a). This is 
much less sensitive to the correlated errors that affect the 
Schechter LF parameters. To do so, we define 'giants' as 
galaxies brighter than Mf606w = — 18 and 'dwarfs' as galax- 
ies with -18 < M F 606W < -15. The DGR - the ratio of the 
numbers of so defined dwarfs to giants - was calculated for 
each dwarf galaxy by finding the distance to its 10th near- 
est giant galaxy and counting the number of dwarf galaxies 
within the same region. The area of the circle was calculated 
numerically to take into account gaps in the mosaic (Figure 
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Figure 10. Left panel: Galaxy luminosity distributions (LDs) 
for A2218, split on the basis of clustercentric radius. The LD 
for galaxies in the inner 'core' region of the cluster (100 < r < 
300 kpc) is represented by circles; those that are filled have had no 
correction for the effects of gravitational lensing, whereas the open 
circles include this correction. The LD for galaxies in the outer 
'halo' region of the cluster (300 < r < 750 kpc) is represented 
by squares, differentiated in the same way to show the effects of 
lensing corrections. Shechter function fits to the uncorrected data 
are shown as the thick and thin solid lines for the inner and outer 
regions, respectively. Right panel: The la, 2a and 3cr error ellipses 
for the Schechter function parameters, A^pg g W (horizontal axis) 
and a (vertical axis), derived from the fits shown in the left-hand 
panel; inner region (thick lines), outer region (thin lines). 



1). The local galaxy density (Dressier 1980) was then simply 
taken as 10 divided by this area. In order to reduce errors, 
the data were then binned according to the local giant galaxy 
density measured in the vicinity of each dwarf galaxy. Field 
subtraction was performed using the smoothed counts from 
Figure 5. 

Figure 11 shows the DGR plotted as a function of lo- 
cal giant galaxy density. Our results show clearly that the 
DGR in A2218 decreases smoothly and monotonically with 
increasing local giant galaxy density, consistent with the 
dwarf population-density relation proposed by Phillipps et 
al. (1998). Indeed our data define the relation much more 
cleanly and over a large range in galaxy density than the 
data from Smith et al. (1997) and Driver et al. (1998a) 
(which we also plot in Figure 11 for a qualitative compar- 
ison) upon which Phillipps et al's proposed relation was 
based. That these latter data points all appear to sit lower 
on the diagram with respect to our A2218 data is most likely 
attributable to differences in niters, magnitude depths and 
binning in galaxy density between the different studies. 

Our data are of sufficient depth to allow us to also 
consider a fainter sample of 'ultra-dwarf galaxies, with 
— 13.5 < Mf606W < —15. The DGR values for this popu- 
lation are also plotted in Figure 11 and appear to follow 
the same trend with density as the brighter sample, indicat- 
ing commonality between the ultra-dwarf and dwarf popu- 
lations. 



3.3 The radial profile of galaxies in A2218 

In order to determine whether or not the decreasing DGR 
is due only to the increasing giant galaxy density toward 
the cluster core or whether there is a corresponding de- 
crease in the number of dwarfs for higher giant galaxy den- 
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Figure 11. The variation of the dwarf-to-giant galaxy ratio 
(DGR) with local giant galaxy density. Filled circles represent the 
'intermediate' dwarfs (—18 < Mpgoew < —15) in A2218; starred 
points represent the 'ultra-dwarfs' (—15 < Mf606W < —13.5) 
in A2218. The error bars include Poisson contributions only and 
lensing effects arc ignored. The open triangles show the data for 
A2554 from Smith et al. (1997). The filled squares represent the 
central regions of the clusters from Driver et al. (1998a) and the 
open squares the outer regions. 



sities, we have derived radial density profiles for these dif- 
ferent sub-populations within A2218. These are plotted in 
Figure 12, with the data for the giants (Mpeoew < —18) 
shown in the top panel, those for the 'intermediate' dwarfs 
(—18 < Mf606w < —15) shown in the middle panel, and 
those for the 'ultra-dwarfs' (—15 < Mpeoew < —13.5) shown 
in the bottom panel. Again the data with and without correc- 
tions for lensing are plotted, and it can be that the difference 
between them is negligible. Accordingly, in fitting the data 
with a power-law (a(r) oc r a ; solid lines in Fig. 12), we do 
so to the uncorrected data only. Note that we also exclude 
the innermost radial point in these fits as well (correspond- 
ing to the inner 100/i -1 kpc where the lensing correction is 
uncertain. 

The giant galaxy projected density is seen to drop off 
with a power-law index a — —1.13 ± 0.15, consistent with 
that of a projected isothermal sphere (oc a -1 ). We also note 
that beyond a radius of ~ 1 arcmin (~ 128 kpc), it traces 
very well the projected X-ray gas density profile (dashed 
line) observed for A2218 (Machacek et al. 2002), which has 
a core radius of r c = 66.4" and f3 — 0.705. 

For the intermediate luminosity dwarfs, the distribution 
is flatter than that of the giants and is best fit by a power- 
law with index a = — 0.63T0.09. In contrast, the ultra-dwarf 
galaxy surface density shows a deficiency in the centre of 
the cluster, and then is approximately flat at r > 0.6 arcmin 
(~ 200ftT 1 kpc). A power- law fit to the data has a index 
of a = —0.03 ± 0.13. The most likely explanation for this 
behavior is that the fainter galaxies have a more extended 
distribution, and hence larger scale radius, than the giants, 
as also seen for the bright and faint dEs in the Coma cluster 
(Seeker et al. 1997). 

There are two competing effects which should be consid- 
ered in interpreting Figure 12. Firstly, because of the larger 
galaxy density in the centre of the cluster there is a dimin- 
ishing area effect where the field of view in which dwarf 
galaxies can be detected is reduced due to the presence of 
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Figure 12. Radial surface density distributions for the giant, 
dwarf, and ultra-dwarf galaxy populations in A2218: Top panel: 
the giant (Mp606W < —18) galaxies, with the uncorrected data 
shown as filled circles, and the lcnsing-corrcctcd data shown as 
open circles. The projected X-ray gas density profile is shown as 
the dashed line. Middle panel: the intermediate dwarf (—18 < 
A^F606W < —15) galaxies, with the uncorrected data shown 
as filled squares and the lensing-corrected data shown as open 
squares. Bottom panel: the ultra-dwarf (—15 < Afpgoew < —13.5) 
galaxies, with the uncorrected data shown as filled trianqles and 
the lensing-corrected data shown as open trianqles. A power-law 
fit to the uncorrected data points is represented in all three panels 
by the solid line. 

brighter galaxies (see Driver et al. 1998b) . To determine the 
extent of this effect, we sum together the areas assigned to 
detected objects by SExtractor in both the very central 
region of the cluster and in the lowest density regions. We 
find in the central region the total area associated with ob- 
jects by SEXTRACTOR is 8.2% of the total area, while in the 
outskirts only 4.2% of the area is occupied by objects. Hence 
if the data in Figure 12 were to be corrected for this effect, it 
would amount to only a 4% differential adjustment between 
the innermost (highest density) and outermost (lowest den- 
sity) bins. This would have a negligible effect in terms of the 
overall trends seen in the bottom two panels. 

Secondly, there is a radial dependence on the volume of 
cluster sampled per unit area projected on the sky. Assuming 
the cluster fills a sphere of radius R, then the volume per 
unit area projected onto a circle of radius a centered on the 
cluster is : 

Volume ^ 4 {R s _ {R * _ a *f 2) (15) 
Area 3a 2 v v ; ; v ' 

Hence the volume per unit area decreases with radius. If we 
assume a physical radius for A2218, then we can calculate 
the corrections in number needed in each of the radial bins in 
Figure 12. For a spherical cluster with a 1.5 /i" 1 Mpc radius 
(approximately the Abell radius), a factor of 1.15 increase in 
number is required in the outermost annuli (with respect to 
the innermost annuli) to convert from number per unit area 
to number per unit volume. This more than compensates for 



the correction (in the opposite sense) due to the diminishing 
area effect. So it is likely that the differences in profile shape 
for the giant, dwarf, and ultra-dwarf galaxy populations are 
as large, if not larger than that seen in Figure 12, and the 
flattening of profile shape with decreasing luminosity is a 
real effect, when both these effects are taken into account. 



4 SUMMARY AND DISCUSSION 

We have exploited deep, wide-field HST imagery of the well 
known rich cluster A2218 at z = 0.18, to conduct a census 
of its galaxy population over a ~ 10 magnitude range in 
luminosity (—23 < Mf606W < —13), and thereby explore 
how its constituent 'giant' and 'dwarf galaxy populations 
vary in number as a function of location and environment 
within the cluster. 

Taking the complete ensemble of A2218 galaxies 
that are identified, statistically, within the entire central 
1.34 Mpc 2 hT 2 field studied and over the full range in lumi- 
nosity and surface brightness (hfsogw g 24.7 mag arcsec -2 ) 
probed, we find their luminosity distribution to be well fit- 
ted by a Schechter function, with parameters Mp 606W = 
-20.52 ± 0.27 and a = -1.38 ± 0.05. In terms of the faint- 
end slope, at least, this is in reasonable agreement with the 
LFs measured for other rich clusters, both at similar and 
more nearby redshifts. 

Importantly, however, our study has also revealed that 
such 'global' measures of the LF represent the superposition 
of a number of important underlying dependencies on clus- 
tercentric radius and, more fundamentally, the local galaxy 
density. The most conspicuous of these is a change in the 
shape of A2218's LF with clustercentric radius, in particu- 
lar its faint end having a steeper slope at larger radii. This 
indicates that the relative balance between the numbers of 
dwarf and giant galaxies varies with environment in the clus- 
ter, and we have quantified this by measuring the ratio of 
dwarf to giant galaxies (DGR) as a function of local pro- 
jected galaxy density. We find the DGR in A2218 to decrease 
monotonically with increasing local density, dropping by an 
order of magnitude (from ~ 10 to ~ 1) in going from the low- 
est (~ 100 gals Mpc" 2 ) to the highest (~ 1000 gals Mpc" 2 ) 
densities probed by our observations. This confirms the exis- 
tence of a "dwarf galaxy-density relation" as first suggested 
by Phillipps et al. (1998). Furthermore, we have shown that 
this relation arises from galaxies of different luminosity hav- 
ing quite different radial surface density profiles. The bright- 
est, giant galaxies have the steepest profile, consistent with 
an isothermal sphere and that of the hot X-ray gas. The 
fainter, dwarf galaxies have much shallower profiles, indi- 
cating that they are a much more extended population with 
a much larger scale radius. Indeed the 'ultra-dwarf galaxy 
population in A2218 is found to have a radial profile which 
is essentially flat, indicating a constant surface density at 
least over the central 750 kpc studied here. 

This quite clear 'segregation' between galaxies in these 
different luminosity regimes could be a result of primordial 
conditions, evolution, or both. As has been previously dis- 
cussed by Phillipps et al. (1998), the numerical CDM models 
of Kauffmann et al. (1997) point to the dwarf population- 
density relation being a natural consequence of the hierar- 
chical clustering process. In this picture the dwarf galaxies, 
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which form via the gravitational collapse of small amplitude 
density peaks in the primordial universe, are predicted to be 
more numerous and less clustered than their brighter coun- 
terparts. 

However, Moore et al. (1998) have also been able to 
explain the dwarf population-density relation as the result 
of galaxy 'harassment', which operates more effectively in 
denser regions. As the cluster tidal field becomes increas- 
ingly important in the cluster centre, the smaller spheroidal 
galaxies will be destroyed in the innermost part of the cluster 
and global tides will disintegrate the lowest surface bright- 
ness objects into the diffuse stellar background. In the Coma 
cluster, for example, Biviano et al. (1996) find that the clus- 
ter structure is better traced by the faint galaxy population, 
which forms a single smooth structure, with the brighter 
galaxies being located in sub-clusters. They interpreted this 
as evidence for the ongoing accretion of groups onto the 
main body of the cluster. In this context, it is interesting to 
note that Bernstein et al. (1995) find a significant flattening 
of the density profile of faint galaxies in the inner 100 kpc of 
the Coma cluster. 

The deficiency of dwarf galaxies in the centre of Abell 
2218 may also be related to the presence of a cD galaxy. 
cD galaxies probably form suddenly during the collapse and 
virialization of compact groups or poor clusters which then 
merge with other groups to form a rich cluster (Merritt 
1985). Lopez-Cruz et al. (1997) suggest that the flatter faint- 
end slope of the galaxy LF which they find in rich clusters, 
in particular clusters that contain cD galaxies, results from 
the disruption of large number of dwarf galaxies early on in 
the clusters evolution. A similar effect is claimed by Bark- 
house et al. (2002). The stars from the disrupted galaxies 
are redistributed throughout the cluster, generating the cD 
halo. This is consistent with a lack of dwarf galaxies in the 
core of Abell 2218. 

Ultimately, a single cluster is not sufficient to both es- 
tablish the trends we discussed and analyse their dependence 
on cluster properties. A larger sample of clusters with deep, 
wide-field imaging is needed to make progress on these im- 
portant issues. This is a project we are currently embarking 
upon, the results from which will be reported in future pa- 
pers. 
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